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Abstract. We describe large scale structure at high redshift in terms of the Cos- 
mic Web picture for {S, A, O, H}CT)M models: how galactic-scale "peak-patches", 
filaments and membranes create an interconnected intergalactic medium. The ideas 
are applied to our Lya forest simulations of "shear-field patches". We discuss sim- 
ulation method and design, resolution dependence, the statistical combination of 
patches, UV flux scaling, and whether filtered Zel'dovich maps are useful. The re- 
sponse to changes in power spectrum shape and amplitude, and in cosmological 
parameters, is described. We also show f2(,h 2 derived from UV rescaling is overesti- 
mated if the resolution is not adequate. 

1 The Cosmic Web and Peak-Patches at High Redshift 

The Cosmic Web picture |lj provides a powerful language for understanding 
the structure and evolution of Lyman Alpha absorption systems Q . It predicts 
the basic structural components of the IGM as a function of scale, epoch and 
cosmology. For density contrasts 6 >100, the rare-events at z ~ 3 are massive 
galaxies, observed as damped Lya absorbers, while the typical collapsed ob- 
jects are dwarf galaxies responsible for Lyman Limit and metal line systems. 
The medium is visually dominated by 5 ~ 5 — 10 filaments, bridging massive 
galaxy peaks, with smaller scale filaments within the larger scale ones bridging 
smaller dwarf galaxies contributing the most to the Njji ^ 10 14 ' 5 Lya forest. 

A-body calculations cannot currently cover large enough volumes of space 
with sufficient resolution to simulate high redshift galaxy catalogues (with 
clustering) to compare with e.g., the Steidel et al. structure at z ~ 3.1 ||. 
Only intermediate and small scale structures are accessible to direct gasdy- 
namical simulations of the Lya forest. A peak-patch catalogue with ultralong 
waves included can treat such large volumes quickly and accurately Q] . Galaxy 
halos are identified with regions of space determined to have collapsed using 
ellipsoidal internal dynamics, with external tidal fields playing a large role. 
Bulk properties like mass, binding energy, velocity and shear for the patches 
follow, and can be used with single-patch hydrodynamics or phenomenology 
to predict internal gas profiles. Field realizations constrained to have inter- 
esting multiple peak/void structures can be used as initial conditions for high 
resolution numerical simulations, e.g. the strong filament of galaxies in 

In Fig. [j], we show the (impressively grand) large scale clustering in galaxies 
at various velocity dispersion cuts, in comoving space and in redshift space. 



Motivated by 0, we tiled regions from z=2.8 to 3.5 encompassing 18' x 18' 
with 128 3 , 40 Mpc boxes. The size was chosen to resolve "dwarflet" peak- 
patches with binding energy vbe ~ 30 km s _1 . Depending upon cosmology 
(Table Q), 15 to 22 boxes were needed. (Tiling for open (hyperbolic) universes 
is an interesting exercise.) Optimal wavenumber sampling was used (§ ||), with 
phase coherent ultralong and short waves consistently joining box to box. 



Table 1: The cosmologies studied are normalized to fit the COBE data, except for 
SCDM. The {S, A, O, //}CDM hydro runs (first parameters given) had Gaussian- 
smoothed ct p (0.5 Mpc) =1.05 at z=3, and O fl h 2 =0.0125. 



Cosmology 


OS 


H , Age 




Notes 


SCDM 


0.67,0.44,1.0 


50, 13 Gyr 


1.0 




ACDM 


0.93, 0.91 


70, 13 Gyr 


0.335 


ja A =0.665,n s = 0.94, n s = 1 


OCDM 


0.81, 0.91 


70, 11 Gyr 


0.37 


n s = 1 


HCDM 


0.81 


50, 13 Gyr 


0.8 


Qmv = 0.2, 2 degenerate species 



Our "shear-patch" gas simulations sample the Fig. [TJ space only over small 
high resolution cores 5 Mpc across (encompassing a fixed cosmology-independent 
baryonic mass of M& = 2 x 10 11 Mg), but our images of Njji,Tg as etc. (see 
p^[ ) reveal the cosmic web: In the most probable shear-patches, the typi- 
cal IGM at z ~ 3 is dominated by filamentary structures fixed in comoving 
space with large dwarf galaxies (M ~ 10 10 ' 5 Mq, vbe ~ 90km s _1 , Nhi ~ 
10 17 ~ 19 cm~ 2 ) forming at the junctures and smaller dwarfs [M ~ 10 9,2 -M Q , 
vbe ~ 30 km s _1 , Nhi ~ 10 15_17 cm -2 ) littering the limbs. In less probable 
protogalactic shear-patches, the filaments merge, as part of larger scale fila- 
ments. In protovoid shear-patches, the webbing is less pronounced. ACDM 
and OCDM have fatter ribbon-like filaments with fewer dwarflets compared 
to SCDM. 



2 Lya Forest Simulations of "Shear-field Patches" 

Optimal fc-space sampling: The flatness of the power spectrum, dcr^/dlnk, 
for wavelengths A — 2irk~ 1 < 6 Mpc means that small periodic boxes must 
exclude substantial large scale power, and can probe only a restricted set of 
shearing patches. We use a hybrid method for applying waves to the initial 
conditions, with special attention to the correct inclusion of long waves Q: an 
FFT is used for short waves, but slow direct transforms are used with power 
law and log sampling for low k, with switching defined by a fc-space volume 
minimization criterion. 

Importance Sampling: Since a single small box is a poor sample of the 
universe, we use a number of shearing patch simulations for each cosmology to 
be tested, each parameterized by a smoothed overdensity (i/ct p (0.5 Mpc)) and 
anisotropic shear tensor. The Gaussian filter i?/=0.5 Mpc subtends a galaxy 
mass, 3 x 10 11 r2„ r (2h) 2 M©. We have found that using a relatively small 
sample of v : { — 1.4, —0.7, 0, +0.7, +1.4}, with the anisotropic shears taken at 
their mean values given is, gives adequate statistical convergence on quantities 
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Figure 1: The redshift-space (left) and comoving-space (right) distributions of galac- 
tic peak-patches. Top left is 18' wide and 9' thick, with a vbe > 290 km s _1 cut, 
roughly corresponding to Steidel's damped Lyce systems ||. The 4 cosmologies 
shown, {S, A, O, //}CDM, all reveal large scale features, though more pronounced in 
OCDM and ACDM. Top right shows the spatial distribution, this time for bright 
galaxies with a 200 km s _1 and 3 x 10 10 Mq baryonic mass cut, for a region now 
4.5' thick. Bottom left and right show a 90 km s _1 cut, for a 1' thick slice. The 
filaments bridging these dwarf galaxies define the lower Nhi Lya forest, while ha- 
los with vbe ~ 30 km s _1 dominate the Nhi ~ 10 15 regime. The forest hydro 
simulations cover only 1/8 of the 40Mpc at high resolution. This figure therefore 
emphasizes large scale waves must be included in our small scale "shear-patch" sim- 
ulations, whose size is governed by our need to resolve 1 kpc structure in galaxies at 
z=3. The patches then typically include many 30 km s" 1 halos but only a handful 
of 90 km s _1 ones. The currently largest periodic simulations with this resolution 
are not much bigger, but miss the long waves. We can use these peak-patch/cosmic 
web simulations to compare predictions for different cosmologies with the large scale 
structure probed by multiple quasar line-of-sight data and long range velocity space 
correlations in quasar spectra as well as emerging high- z catalogues M] . 



0.035 - 



0.030 



0.015 







7 Cen ACDMo 


■ 


■ 

~ Low Res 


SCDM - 


SCDM^ 


■ 


" Croft SCDM^/ 


cy„=0.7 




o 8 =0.4 : 








cr 8 =o.7 ■ 




H & M 



1.5 2.0 2.5 



3.0 3.5 
Redshift Z 



4.0 4.5 5.0 



Figure 2: Qt h 2 derived from UV flux rescaling for the 4 cosmological models at z=3, 
and as a function of redshift for SCDM for the 3 <r p (0.5 Mpc) values we tested. A fixed 
J_2i=0.5, typical of z~3 estimates, is used, except for the evolving H&M (SCDM) 
curve using the Hardt and Madau UV flux Q. Low Res SCDM is a simulation 
with a factor of 3 lower resolution than our standard, which gives incorrect opacity 
estimates. "Cen" and "Croft" refer to simulations discussed in EJ and the text. 
Raising as (i.e., <7p(0.5 Mpc)) results in a lower mean opacity by emptying high 
cross section filaments into more compact objects. J-21 is the flux at the hydrogen 
edge in 10 -21 cgs units. 
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Figure 3: Column density distribution at z=3, with a iV^}' 46 power-law removed, 
reveals structural details, including the "break" at Nhi ~ 10 14 ' 5 cm -2 . Results below 
10 16 cm -2 fit the observations well for all 4 cosmologies. This is partly due to the 
forgiving nature of opacity rescaling that prevents differentiation among cosmological 
models based on this measure, even when spatial images of the simulations in Nhi 
look very different. The trend with redshift [[lo| is also reproduced. 



such as the frequency distribution of Nhi , provided the simulations are given 
their proper relative weighting. The weighting of ^-bin i is oc J u ' +1/2 eT v l 2 dv 
in initial condition space, but in the dynamically evolved space it must be 
multiplied by the angle on the sky subtended by the region and the length in 
redshift space. The negative v (void) cases expand and thus contribute more, 
reducing the frequency of lines in the combined sample relative to what an 
average density simulation would give. Overdense regions are more compact 
and contribute large numbers of high column lines to the statistics (e.g., the 
v=\A protogalactic patches form collapsed massive galaxies just below z ~ 2). 
SPH-TreeP 3 M. The simulations reported here used a new code, combining 
an FFT for the P-M [|| with a tree-like P-P, which is twice faster than the 
SPH-P 3 MG code used in our earlier work and with more accurate forces. 
Our high resolution 5.0 Mpc diameter patches had 50 3 /2 gas and 50 3 /2 dark 
matter particles (initial grid spacing of 100 comoving kpc), surrounded by gas 
and dark particles with 8 times the mass to 8.0 Mpc, in turn surrounded by 
"tidal" particles with 64 times the mass to 12.8 Mpc. A linearly-evolved self- 
consistent external shear field was included to treat ultralong waves. SPH is 
adaptive: our best resolution is set by the gravitational softening of 1 kpc. 
Numerical convergence is fundamental but has been largely neglected in 
previous studies. Simulations with the same initial conditions, but 100, 50 
and 25 comoving kpc initial spacing, show our standard 100 kpc computa- 
tions are close to convergence in the Nhi images and in the column density 
distribution /(Nhi), but of course not in small scale dark matter clumping. 
However, tests halving the initial resolution to 200 kpc initial spacing show 
a significantly lowered mean opacity and fail to produce small dwarf galaxies 
(removing objects below ~ 10 10 M Q ). 

3 Lya Forest Response to Varying Physical Parameters 

We have varied cosmology, UV flux history and power spectrum shape and 
amplitude Jl3| while keeping the random number choices for the initial wave 
phases and the simulation size fixed, to facilitate comparisons. 
The flux decrement one-point distribution: P(> Da), Da = 1 e r , 
is determined mainly by the mean opacity and is insensitive to choices of 
cosmology. Fitting the observed P(> Da) well after rescaling is all too easy; 
we have even shown that P(> Da) data from other redshifts can be rescaled 
to fit it. (Although we use the full P(> Da) to rescale, just scaling to < Da > 
gives identical results.) 

UV flux history: We have tested how well the much-used opacity rescaling 
by (fish 2 ) 2 /( J- 2 iff (z)) actually works, by evolving identical patches except 
for the flux history. We find it works quite well in the intercloud filamentary 
regime but can fail in halos (the traditional clouds), where the Jeans Mass 
history is important, and thus can fail for the higher column lines.^j 

We have also shown the truncated Zel'dovich map with a global filter Rf and 



fib and UV flux: Fig. g shows the scaling parameter for the four cosmologies 
for J_2i = 0.5. At z = 3, (ft b h 2 )/(2J_ 2 i) 1/2 ranges from 0.018 for ACDM 
to 0.020 for OCDM. Note that our one-third-resolution results (labeled Low 
Res SCDM) give a significantly higher estimate; the estimate labeled Croft 
[ |l~2| also had about this resolution, partly explaining why it is high, but also a 
colder thermal history. We also note that the higher flux or lower Oj, needed 
to scale ACDM that we see follows from the slower expansion rate at high 
redshift than SCDM has (cf. the Cen result ||). 

Column density distribution: Fig. || shows that with the rescaled UV flux, 
J(Nhi) is well-fit by all 4 cosmologies tested. They are visually different, 
which higher-point statistics (e.g. 2-point (^]) can better reveal. The f(Nni) 
were generated from automated Voigt profile fits to simulated spectra. Each 
curve represents the weighted composite of the 5 v simulations. Our simulated 
results do exhibit an evolutionary trend in the Lya line counts similar to 
the Keck observations of jnl, a consequence here of the shift in the knee 
in the distribution apparent at 10 14 ' 5 cm~ 2 , which marks the transition from 
filamentary IGM material to the halos of dwarf galaxies, 
us Variation: For SCDM, we varied the overall power spectrum amplitude: 
cr p (0.5 Mpc, z = 3)=0.7,1.58 as well as 1.05. Although dramatically different 
in visual appearance at the same redshift, with UV rescaling the f(Nni) dis- 
tributions look rather similar at low Nhi- The required scaling is larger than 
for the cosmology variations with fixed <r p (0.5) of Fig. 0. The median in the 
line width, b, is ~ 2 km s _1 higher for 1.58, and ~ 2 km s _1 lower for 0.7, cf. 
the 1.05 case (with median « 30 km s" 1 , as for {A, O, H}CDM). 
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<jp(Rf) < 1, gives images which differ radically from our simulations. Allowing Rf to 
vary reproduces }{Njjj) for low Njjj, a testament to the forgiving nature of the f(Njjj) 
distribution with UV flux rescaling {cf. 



